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Abstract. The angular distributions for elastic scattering and breakup of halo nuclei are 
analysed using a near-side/far-side decomposition within the framework of the dynamical eikonal 
approximation. This analysis is performed for n Be impinging on Pb at 69ylMeV. These 
distributions exhibit very similar features. In particular they are both near-side dominated, 
as expected from Coulomb-dominated reactions. The general shape of these distributions is 
sensitive mostly to the projectile-target interactions, but is also affected by the extension of the 
halo. This suggests that the link between elastic scattering and a possible loss of flux towards 
the breakup channel is not obvious. 



1. Introduction 

The development of radioactive-ion beams in the mid-80s has enabled the exploration of the 
nuclear landscape away from stability. This technical breakthrough led to the discovery of 
exotic nuclear structures such as haloes [HE]- Halo nuclei are light neutron-rich nuclei, which 
exhibit a matter radius significantly larger than their isobars. This large size is qualitatively 
understood as resulting from their small binding energy for one or two neutrons [3]: Due to 
their loose binding, these valence neutrons can tunnel far away from the core of the nucleus 
and exhibit a large probability of presence at a large distance from the other nucleons. Halo 
nuclei have thus a strongly clusterised structure: they can be seen as a core to which one or two 
neutrons are loosely bound. These valence neutrons hence form a sort of diffuse halo around a 
compact core. The best known halo nuclei are n Be and 15 C, with a one-neutron halo, and 6 He 
and 11 Li, with a two-neutron halo. Proton haloes can also develop around proton-rich nuclei, 
such as 8 B or 17 F. 

Since their discovery, halo nuclei have been at the centre of many studies, both experimental 
[U [5] and theoretical [6j \7\. Due to their short lifetime, they cannot be studied with usual 
spectroscopic techniques, and one must resort to indirect methods to deduce information about 
their structure. Reactions are the most used tools to study halo nuclei. In particular, elastic 
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scattering [9] and breakup [TOJ E] convey interesting information about the structure of the 
projectile. 

Recent experimental [9] and theoretical [8] studies of elastic scattering indicate a strong 
coupling between scattering and breakup. On the experimental side, the elastic scattering cross 
section for 11 Be on Zn around the Coulomb barrier is significantly reduced at large angles 
compared to that of non-halo Be isotopes [9j. One explanation of this unexpected reduction 
is the transfer of probability flux from the elastic channel to the breakup channel. On the 
theoretical side, Matsumoto et al have shown that for the elastic scattering of 6 He on Bi at low 
energy CDCC calculations agree with experimental data only if the breakup channel is included 
in the model space. 

To better investigate the interplay between elastic scattering and breakup, we analyse 
theoretically the angular distributions for the elastic scattering and breakup of halo nuclei within 
a near/far decomposition |12t 113] . We choose 11 Be, the archetypical one-neutron halo nucleus, 
as testcase and consider its collision on Pb at 69^4MeV, which corresponds to the conditions 
of the RIKEN experiment [TO]. The calculations are performed within the Dynamical Eikonal 
Approximation (DEA) |1H I14j . which is in excellent agreement with various experimental results. 

After a brief reminder of the DEA and the near /far decomposition, we apply this technique 
to the elastic-scattering cross section (Sec. [3]). We then move to the analysis of the angular 
distribution for breakup (Sec. H|) and show how similar both cross sections are at intermediate 
energies. In Sec. we emphasise the consequences of this analysis for the study of halo nuclei 
and provide the prospects of this work. 

2. Theoretical framework 

2.1. Dynamical eikonal approximation 

Most of the models of reactions involving one-neutron halo nuclei rely on a three-body description 
of the colliding nuclei [UJ [TS]: a two-body projectile P made up of a fragment / loosely bound 
to a core c impinging on a structureless target T. The two-body structure of the projectile is 
described by the phenomenological Hamiltonian 



where r is the c-f relative coordinate, /x c j is the c-f reduced mass, and V c f is a real potential 
adjusted to reproduce the binding energy of the fragment to the core and some of the excited 
states of the projectile. This potential usually exhibits a Woods-Saxon form factor and may 
include a spin-orbit coupling term. 

The interaction between the projectile components c and / and the target T are simulated 
by the optical potentials V c t and VfT, respectively. Within this three-body framework, studying 
reactions involving one-neutron halo nuclei reduces to solve the three-body Schrodinger equation 
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A R + H + V cT (r, R) + V fT (r, R) (r, R) = E tot V(r, R), 
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where R is the P-T relative coordinate, \x is the P-T reduced mass and 
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is the total energy of the system, with E$ the (negative) energy of the projectile ground state 
<pl j m an d UK the initial P-T relative momentum. The quantum numbers Iq, jo and mo 



correspond to the c-f orbital angular momentum, the projectile total angular momentum and 
its projection, respectively. 

To describe a reaction in which the halo nucleus P impinges on the target T, Eq. ([2]) is solved 
with the initial condition 

^(r,b,Z) — > e iKZ ct> lojomo , (4) 

Z— i>— oo 

where the dependence upon the transverse b and longitudinal Z components of R is made 
explicit. Equation ([2]) must be solved for each value of b and of uiq. 

At sufficiently high incident energy, the eikonal approximation can be performed to ease the 
resolution of Eq. ([2]) . That approximation consists in assuming that most of the rapid variation 
of in the P-T relative coordinate Ft is included in the plane wave e iKZ , i.e. that the three- 
body wave function is well approximated by that plane wave times a function ^ that does not 
vary much with R: 

*(r,b,Z) = e tKZ $(r,b,Z). (5) 
Including the eikonal ansatz ([5]) within Eq. ([2]) leads to 
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V(r,R)=E tot *(r,R). (6) 



Since varies smoothly with R, its second-order derivative A^\& can be neglected in front of its 
first-oder derivative Kd/dZ^f. Then, considering the energy conservation ([3|), the three-body 
Schrodinger equation ([6]) reduces to the DEA equation [TTJ [T3] 

ihv^(r,R) = [H -E + V cT (r,R) + V fT (r,R)]^(r,R), (7) 

with v = fiKj [i the initial P-T relative velocity. This equation is mathematically equivalent to a 
time-dependent Schrodinger equation with straight-line trajectories posing Z = vt. It can thus 
been solved using appropriate algorithms, such as the one described in Ref. [16] . However, since 
the DEA does not assume any semiclassical treatment of the P-T relative motion, b and Z are 
quantal variables. This implies that the DEA includes quantal interferences such as between 
different trajectories, which are missing in time-dependent models [13] • The DEA therefore 
significantly improves these models. 

The DEA differs also from what is usually called the eikonal approximation [T5] . In its usual 
form, the eikonal approximation includes a subsequent adiabatic approximation to Eq. ([7]) in 
which the excitation energy of the projectile is neglected, i.e. H$ — Eq ~ 0. In that case, the 
solution of Eq. is approximated by the eikonal form factor 

§&° ) (r,6,Z oo) = e^ r ' b Ui 03omo (r), (8) 
where the eikonal phase reads 

X(r, b) = -^~ J™[V cT (r, b, Z) + V fT (r, b, Z)]dZ. (9) 



The DEA thus improves the usual eikonal approximation by including dynamical effects that 
are otherwise neglected. These effects may be very significant such as in Coulomb breakup, for 
which the usual eikonal approximation diverges [14] . 



The DEA has been used successfully to describe elastic scattering and breakup of one-neutron 
halo nuclei on both light and heavy targets |14| . This approximation has also provided a reliable 
description of the Coulomb breakup of the one-proton halo nucleus 8 B [17] . More recently, a 
comparison of various reaction models has shown that the DEA is in excellent agreement with 
CDCC at intermediate energies [18] while being much less demanding on a computational point 
of view. The DEA is thus the most efficient model to study reactions involving one-neutron halo 
nuclei at intermediate energies. Moreover as it describes simultaneously both elastic scattering 
and breakup, the DEA is ideal for the present study. 



2.2. Angular distributions and their near/far decompositions 

Within the DEA, the angular distribution for elastic scattering, i.e. the elastic-scattering cross 
section, reads [2] 
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where is a Bessel function [19] . q = 2K sin 6 /2 is the transferred momentum and 
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To have a better insight into the reaction mechanism that takes place during the scattering 
of the projectile by the target, we perform a near/far decomposition of the elastic-scattering 
cross section (|10p |12[ 113]. The idea behind this decomposition is to express the Bessel function 
as the sum of two Hankel functions 
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The elastic-scattering cross section can then be expressed as the sum of two terms obtained by 
substituting J M by either H^' /2 or /2 in Eq. (|10|) . The former is called the Far side (F) of 
the angular distribution ([10]) , while the latter is its Near side (N): 
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Since these two terms add coherently to form the elastic-scattering cross section, they may 
interfere when they reach similar magnitude, which explains some of the oscillatory patterns 
observed in angular distributions [13] . 

The physics behind this decomposition can be understood from the asymptotic behaviour of 
the Hankel functions: 
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Since q K9, the N side corresponds to the positive deflection i.e. repulsive forces (see Fig. [I]). 
On the contrary the F side carries information about negative deflection, i.e. attractive forces. 

The angular distribution for the breakup of the projectile can also be computed within 
the DEA [33]. It corresponds to the breakup cross section expressed as a function of the 
scattering angle Q = (9, (p) of the c-f centre of mass after dissociation at a c-f relative energy 
E = h 2 k 2 /2/j, c f. It reads 
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Figure 1. Schematic illustration of the Near and Far sides of the angular distribution, 
where 

Sfi$(«0 = (0ki jm \* imo \b,Z -> oo)>, (16) 

with 4>kijm the continuum wave function describing the broken up projectile in partial wave Ijm. 

To study the breakup process, we extend the N/F analysis to the angular distribution (fT5|) . As 
for the elastic scattering, the Bessel function is decomposed into the sum of two Hankel functions 
([12p . which provide both N and F sides of the breakup cross section with the same interpretation, 
i.e. the contribution to breakup of the repulsive and attractive forces, respectively. 

3. Elastic scattering 

3.1. n Be on Pb at 69 AMeV 

As a first step in our analysis, we study the elastic scattering of Be on Pb at 69j4MeV. 
As mentioned earlier, 11 Be is the archetypical one-neutron halo nucleus. In our analysis, it 
is described as a 10 Be core in its + ground state to which one-neutron is loosely bound. We 
choose for the 10 Be-n interaction the potential developed in Ref. [20], which reproduces the l/2 + 
ground state in the lsl/2 partial wave at 504 keV below the one-neutron separation threshold. 
This potential also reproduces the l/2~ bound excited state in the Opl/2 orbital and the 5/2 + 
resonance in the d5/2 partial wave. We use the numerical parameters and the potentials V c t 
and Vfr detailed in Ref. [2]. The numerical technique used to solve the DEA equation ([7]) is 
explained in Ref. [TB] 

The DEA elastic-scattering cross section is plotted in Fig. as a ratio to Rutherford [2T]. It 
presents a usual shape with a Coulomb rainbow at about 2° followed by an exponential drop. At 
larger angles, the angular distribution presents significant oscillations. The N/F decomposition 
shows that at forward angles the process is fully dominated by the N side, as expected for 
a (repulsive) Coulomb-dominated reaction [13]. Note that the forward-angles oscillations (i.e. 
below 2°) are observed in both the total cross section and its N side. The N/F interferences 
therefore cannot explain this feature of the angular distribution. At larger angles, i.e. around 
8°, the N and F sides cross, explaining the oscillatory pattern of the total cross section. This 
shows in particular that attractive nuclear forces affect the elastic scattering mostly at large 
angles, as is expected from semiclassical models. 

The whole interpretation of the N/F decomposition is based on the asymptotic behaviour of 
the Hankel functions (|14p . To validate this interpretation, we repeat the calculation of the cross 
section, its N and F sides using the asymptotics of the Bessel and Hankel functions. The angular 
distributions obtained in this way (dotted lines) are in excellent agreement with the exact ones, 
confirming our analysis of this N/F decomposition. 
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Figure 2. N/F analysis of the elastic scattering of Be on Pb at 69^4MeV [21 
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Figure 3. Sensitivity of the elastic-scattering cross section to the P-T interaction (left) and 
the extension of the halo (right). 



3.2. Influence of P-T interaction 

To better apprehend the influence of the choice of P-T interaction on the elastic scattering, 
we repeat the calculation with the sole Coulomb term of the nuclear optical potential (i.e. a 
point-sphere potential, P-S) and a purely point-point Coulomb interaction (P-P). The dominant 
N sides of the corresponding elastic-scattering cross sections are plotted in Fig. [3] (left). This 
change of potentials causes dramatic changes in the angular distribution. It mostly modifies 
the Coulomb rainbow. As the full optical potentials, the P-S interaction leads to a Coulomb 
rainbow, but its location is shifted from 2° to about 4°. On the contrary, no Coulomb rainbow 
is observed with the P-P interaction. This confirms that the features of the elastic-scattering 
cross section strongly depends on the choice of the optical potentials V c t and VfT- 

Interestingly, the change in the elastic scattering cross section cannot be simply related 
to a transfer of flux towards the breakup channel, as postulated in Ref. [9]. Although the 
elastic scattering increases at large angles from the Coulomb + nuclear potential to the purely 
Coulomb interaction (first with P-S and then even more with P-P), the total breakup cross 
section increases as well, as shown in Table [TJ Since the Coulomb rainbow appears only for the 
P-T potentials that account for the extension of the projectile and the target (i.e. the full optical 
potential or just its Coulomb component P-S), we now analyse the influence of the extension of 
the halo on these distributions. 
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P-P 


P-S 


C.+N. 


\Eq\ (MeV) 
o"bu (b) 


0.5 
2.58 


0.5 
2.10 


0.5 
1.70 


0.05 
23.57 


5 

0.07 



Table 1. Total breakup cross sections corresponding to the calculations shown in Sees. I3T21 and 
EO ||2T]. 
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Figure 4. N/F decomposition of the breakup angular distribution for n Be on Pb at 69AMeV. 
The 10 Be-n relative energy is E = 0.5 MeV. 



3. 3. Influence of the size of the halo 

To study the sensitivity of our results to the size of the halo, we repeat the DEA calculations 
adjusting the 10 Be-n potential to increase (reduce) the neutron separation energy \Eq\ of the 
11 Be- like projectile in order to shrink (resp. expand) its halo. The N side of the elastic-scattering 
cross section is plotted as a ratio to Rutherford in Fig. [3] (right). 

The slope of the exponential drop after the Coulomb rainbow is sensitive to Eq: Reducing 
the one-neutron separation energy of the projectile, i.e. expanding its halo, slightly reduces the 
elastic-scattering cross section. This effect could therefore be used to get information about the 
extension of the halo. However, this dependence remains small in comparison to the influence 
of the P-T potential [see Fig. [3] (left)] . There is thus little hope that observing the sole elastic- 
scattering cross section could provide unbiased information about the extension of the halo. 

Since reducing \Eq \ increases the breakup cross section (see Table [1]), we could believe that the 
transfer of flux to the breakup channel explains the change in the elastic-scattering cross section. 
However, since this increase is much more significant than the drop in the elastic-scattering cross 
section, our analysis confirms that there is no direct link between both effects, as suggested in 
Sec. ET21 

4. Breakup of n Be on Pb at 69,4MeV 

To better comprehend the link between angular distributions for elastic scattering and breakup, 
we perform the same analysis as in Sec. [3] for the breakup cross section (|15p . For the breakup of 
n Be on Pb at 69AMeV, the angular distribution (solid line) and its N (short-dashed line) and 
F (long-dashed line) sides are plotted as a function of the scattering angle 6 of the 10 Be-n centre 
of mass after dissociation |21j (see Fig. 2]). The 10 Be-n relative energy is E = 0.5 MeV. 

The features of the breakup angular distribution are very similar to those of the elastic- 
scattering cross section. First, the full calculation exhibits small oscillations at forward angles 
before an exponential drop starting at 2°, which is reminiscent of the Coulomb rainbow observed 
in the elastic-scattering cross section (see Fig. [2]). Second, at forward angles, the breakup is 
dominated by its N side, just as in the elastic scattering. Finally, at larger angles, the total 
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Figure 5. Sensitivity of the angular distribution for breakup to the P-T interaction (left) and 
the extension of the halo (right). 



cross section exhibits oscillations that are explained by interferences between the N and F side, 
which cross at about 8°. 

Fig. [5] shows the sensitivity of the breakup cross section to the P-T potential (left) and to 
the binding energy of the neutron \Eq\ (right). The similarity between elastic scattering and 
breakup is also observed here. Using the sole Coulomb part of the optical potentials (P-S) shifts 
the start of the exponential drop of the breakup cross section to 4°, as in the elastic-scattering 
one [see Fig. [3] (left)], and using the purely point-point Coulomb interaction (P-P) leads to no 
rainbow-like behaviour. Note that the larger breakup cross section obtained in the P-P case 
(see Table [1]) is explained by this absence of Coulomb rainbow in the angular distribution for 
breakup. The sensitivity of breakup calculations to the neutron separation energy is also similar 
to that observed in the elastic channel. In particular for the slope of the drop after 2°, which 
becomes steeper when |.Eo| is reduced [see Fig. [3] (right)]. This confirms that there is no direct 
link between the drop observed in the elastic-scattering cross section and a possible loss of flux 
towards the breakup channel since both angular distributions vary in the same way when either 
the P-T interaction or the projectile structure are modified. 

As shown in Ref. [22J , these similarities can be semi-quantitatively explained within the Recoil 
Excitation and Breakup (REB) model developed by Johnson et al [23]. Within this adiabatic 
model, the angular distributions elegantly factorise into the product of an elastic-scattering cross 
section for a pointlike projectile and a form factor that accounts for the extension of the halo 
|23} I22j. The fact that the former appears in both factorisations and that it contains most of 
the angular dependence explains the similarity between the angular distributions. That analysis 
also suggests a new observable to study more precisely the halo structure. As observed in 
Fig. [3] (left) and Fig. [5] (left), the angular distributions are very sensitive to the choice of optical 
potentials. Since this sensitivity is very similar in both processes, taking the ratio of two angular 
distributions removes most of the dependence on the choice of the P-T potentials. Such a ratio 
emphasises the nuclear-structure content of the angular distributions [22] . 



5. Conclusion and prospect 

In this work, we analyse theoretically elastic-scattering and breakup reactions of halo nuclei 
through a N/F decomposition of their angular distributions [21] . The calculations are performed 
for 11 Be, the archetypical one-neutron halo nucleus, impinging on Pb at 69^4MeV, which 
corresponds to the experimental conditions of Ref. [TU]. The calculations are performed with 
the DEA, a reliable and accurate reaction model in which elastic scatering and breakup are 
described simultaneously [TTJ E] ■ 

Our analysis shows that at intermediate energy, the angular distribution for breakup is very 



similar to the elastic-scattering cross section: Both present a Coulomb rainbow at the same 
scattering angle 6, they are both N-side dominated at forward angles, they both exhibit similar 
sensitivity to the choice of P-T interaction and to the binding energy of the halo neutron. These 
results suggest that the projectile is scattered by the target in a similar way whether it remains 
bound or it is broken up. This can be semi-quantitatively understood within the REB model 
[221 [23]. 

The present analysis also suggests that there is no obvious link between the drop observed 
in the elastic-scattering cross section and a possible transfer of probability flux towards the 
breakup channel, as postulated in Ref. [9]. Since the work of Di Pietro et al has been performed 
at lower energy (around the Coulomb barrier), our conclusions cannot be directly transposed to 
their study. A similar analysis within the CDCC framework [23] is planned to see how elastic 
scattering and breakup are related to each other at low energy. Moreover, recent progresses 
having been made in the modelling of reactions involving two- neutron halo nuclei [25\ [26] , an 
extension of this work for Borromean systems is also planned. 
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